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I. THE RELATIVITY OF THE QUASIPARTICLE CONTENT OF A QUANTUM FIELD 

t^J- ■ The particle content of a quantum field in flat or curved space-time depends on the motional 

CN | state of the observer, on the way the observer's particle detector couples to the quantum field, and 

on the frequency standard in which the detector carried by the observer measures the quanta to 
be detected. A particular manifestation of this observer dependence is the Unruh effect, which 
consists in the fact that a constantly accelerated detector in the Minkowski vacuum responds as 
if it were placed in a thermal bath with temperature proportional to its acceleration [JO]. The 
deceptively simple fact that the particle content of a quantum field, measured by a suitable detector 
attached to the observer in a supposedly empty Minkowski spacetime, depends on the motional 
state of the observer has eluded direct observation so far. The value of the Unruh temperature 
?Unruh = [fr/ (27r£;BCL)]a = 4 K x a[10 20 g ffi ] , where a is the acceleration of the detector in Minkowski 
space (g® is the gravity acceleration on the surface of the Earth), and cl the speed of light, makes 
it obvious that an observation of the effect is a less than trivial undertaking. Proposals for a 
measurement with ultraintense short pulses of electromagnetic radiation have been put forward in, 
e.g., Refs. 0-0] ) cf. for an extensive discussion and more references the review 0|. Apart from the 
smallness of the effect due to a proportionality constant containing the ratio of Planck constant 
divided by Boltzmann constant times speed of light, the tiny thermal radiation signal needs to 
be detected on top of a (in comparison) huge background signal and in the potential presence of 
qq , dissipative effects. 

In analogue gravity the excitations above a given vacuum, propagating in a generally 

curved spacetime, are quasiparticles like, e.g., phonons, antiferromagnetic magnons, and photons in 
nonlinear dielectric media. An analogue gravity setup of observer-related quantum-field-theoretical 
phenomena has one distinct advantage over real, i.e. Einsteinian, gravity: there exists the preferred 
' reference frame "laboratory" , which allows to state in an absolute manner whether actual dissipa- 

tion has taken place, that is whether the quasiparticles created are due to the choice of the observer 
only, or whether some friction effect creates them, slowing down the motion of the system under 
consideration as observed by the experimentalist in the lab. Two archetypical examples for the 
latter dissipative process are cosmological particle creation due to nonadiabatic effects (squeezing 
of the vacuum) and Hawking radiation from black hole horizons. This sets the purely observer- 
related Unruh effect apart from dissipative effects in an operational, that is measurable manner. 
At the same time, the existence of the absolute reference frame laboratory represents a certain 
limitation of analogue gravity, because it implies the (re-)introduction of an "ether". While an 
ether is to some extent implicit in any quantum field theory, which is ultimately a theory of quasi- 
particles above some vacuum [r|] , the ether in question here introduces a preferred "absolute" time 
coordinate (measured by a clock ticking in the lab frame), a concept alien to relativity. 

It rapidly became apparent that the simulation of effective space-times can be very efficiently 
and with high precision carried out in Bose-Einstein condensates (BECs) , where the relevant quasi- 



particles are phonons 13]. BECs offer the primary advantages of dissipation- free superflow, high 



controllability of the physical parameters involved, and the accessibility of ultralow temperatures 
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[181 ] . Analogues of the dissipative phenomena of cosmological particle creation and Hawking ra- 
diation have been discussed extensively in the literature by now cf., e.g. 
particle creation and 2^-25] for Hawking radiation by sonic black holes 



26 



IS -22 1 for cosmological 



271 ] . As regards the 



experimental realization of such "dumb" holes in BECs, some progress has been made as well 
[ia |. The nondissipative Gibbons-Hawking effect f29|. w hich is the analogue of the Unruh effect 
for a freely falling observer in a de Sitter spacetime [301 ] . will be discussed in what follows in its 
analogue gravity implementation for BECs [3ll. l32j]. For the Unruh effect itself, some aspects of its 
possible realization in BECs, using the same type of detector as first employed in the latter two 
references, have been put forward in 33]. We note that the experimental sensitivity for detecting 
phonons in ultracold dilute quantum gases can reach down to the few or even single-phonon level, 
which is a basic requirement for detecting the subtle effects of quantum radiation phenomena [34| . 
We also note that an alternative setup to detect the analogue Gibbons-Hawking effect has been 
discussed in ion traps 35f] (as well as cosmological particle creation and Hawking radiation, see 
[3^ . I37I]). While the analogy to an expanding universe does not hold there, as there is no real 
(spatial) expansion of the ion array taking place, the advantage of such an ion trap setup is that 
the measurement, because of the spatial stationarity of the central ion locations, be made over a 
comparatively long time. The modulation of the laser coupling electronic states of the ions and 
their vibrational motion (constituting the phononic quasiparticles to be detected), there encodes 
the necessary exponential relation between detector proper time and de Sitter time. 



II. EXPANDING BOSE-EINSTEIN CONDENSATES 



In a seminal paper for the very idea of analogue gravity to be established, it was shown by 
Unruh that the action of the phase fluctuations in a moving inhomogeneous superfluid may be 
written in the form 0], also see 12] (we set % = m = 1, where m is the mass of a superfluid 
constituent particle): 
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Here, v(x,t) is the superfluid background velocity, c(x,t) = y / gpo(x, t) is the velocity of sound, 
where g is a constant describing the interaction between the constituent particles in the superfluid 
(1/g is the compressibility of the fluid), and po(x,t) is the background density. In the second line 
of dl]), the conventional hydrodynamic action is identified with the action of a minimally coupled 
scalar field in an effective space-time which is generally curved [2]. Furthermore, the velocity 
potential of the sound perturbations in the BEC satisfies the canonical commutation relations of 



a relativistic scalar field l4 |. The theory of phononic quasiparticles in the inhomogeneous BEC 



is thus kinematically identical to that of a massless scalar field propagating on the background of 
curved space-time in D + 1 dimensions. We therefore have the exact mapping, on this level of 
kinematics, of the equation of motion for phononic quasiparticles in a nonrelativistic superfluid, 
to quantized massless scalar fields propagating on a curved space-time background with local 
Lorentz invariance. While some dynamical aspects of gravity like backreaction and the cosmological 
constant issue have been studied as well in BECs, see e.g. [H,!^, the actual Einsteinian gravity 
part of the action, proportional to the Ricci curvature scalar, is suppressed by the small gas 
parameter D p = (poa^) 1//2 measuring the diluteness of the system 4(3, 41], where a s is the scattering 
length characterizing atomic collisions in the dilute and ultracold BEC, related to the coupling 



3 



constant by g = 47ra s . Hence Einstein gravity will be subdominant in determining the background 
metric as long as the gas is dilute with respect to the interaction range. 

To model effective curved space-times for quasiparticles, we will consider the evolution of the 
BEC if we change the harmonic trapping, keeping the BEC in a well-defined region of space in the 
lab, with time. For a description of the ex pan sion (or contraction) of a BEC, the so-called scaling 
solution approach is conventionally used 42|. One starts from a cigar-shaped BEC containing 
a large number of constituent particles, i.e., which is in the so-called Thomas-Fermi (TF) limit, 
which in BEC terminology eliminates any gradients of the condensate wave function [43| . According 
to [13], the evolution of a Bose-condensed atom cloud under temporal variation of the trapping 
frequencies co\\(t) and u±(t) (in the axial and radial directions, respectively) can then be described 
by the following solution for the condensate wave function 



TF 



b±y/b 



exp 



bz 2 b_\_v 2 
gp (x = 0, t)dt + i— + i-^— 



(2) 



Here, b± and b are the scaling parameters describing the condensate evolution in the radial (r) and 
axial (z) directions, cf. Fig.[TJ The initial (b = b± = 1) mean-field condensate density is given by 
the usual TF expression 



I^tfI 
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Here, p m is the maximum density (in the center of the cloud) and the squared initial TF radii are 
R?, = 2p/ui 2 and R\ = 2/x/u; 2 . The initial chemical potential fi = p m g. In our cylindrical 3D trap, 
we have for the initial central density 
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The condition that TF be valid implies that fi ^> LJn, lj±. The solution ([2]) of the Gross-Pitaevskh 
mean-field equations becomes exact in this TF limit, independent of the ratio lj^/lj±. However, 
the solution becomes exact also in the limit that — > 0, independent of the validity of the TF 
limit, the system then acquiring an effectively two-dimensional character 44j,|45|. We will see below 
that in the limit of a very long cigar, un/u)±_ — > 0, there is an "adiabatic basis" in which no axial 
excitations are created during the expansion. With respect to that basis there are, in particular, 
no unstable solutions possible, implying the stability of the expanding gas against perturbations. 
According to ([2]), the condensate density evolves as 



Po(r,z,t) 



b'b 



(4) 



and the superfluid velocity 



b ± b 
—re r + -ze 2 
o_i_ o 



(5) 



is the gradient of the condensate phase in Eq. ^j. It increases linearly to the axial and radial 
boundaries of the condensate. 
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FIG. 1. Expansion of a cigar-shaped Bose-Einstein condensate as the sonic analogue of a 1+1D de Sitter 
universe. The violet surfaces represent laser beams creating a tight optical potential well in the center, 
hosting the two- level system Atomic Quantum Dot (red), acting as the detector. The dark surfaces at 
symmetric axial locations zh and — zh represent the sonic horizon. The scaling of the expansion rates along 
axial and perpendicular directions with time, b(t) and b±(t), are derived in the text. 

The excitations in the limit uj^/uj± — > were studied in [iH ]. The description of the modes is 
based on an adiabatic separation for the axial and longitudinal variables of the phase fluctuation 
field: 

<f>(r,z,t)=J2<t>n(r)Xn{z,t), (6) 
n 

where <f>n( r ) is the radial wavefunction characterized by the quantum number n (we consider only 
zero angular momentum modes) . The above ansatz incorporates the fact that for strongly elongated 
traps the dynamics of the condensate motion separates into a fast radial motion and a slow axial 
motion, which are essentially independent. The Xn( z ,t) are the mode functions for travelling wave 
solutions in the z direction (plane waves for a condensate at rest read Xn oc exp[— ie n kt + kz]). The 
radial motion is assumed to be "stiff" such that the radial part is effectively time independent, 
because the radial time scale for adjustment of the density distribution after a perturbation is 
much less than the axial oscillation time scales of interest. The ansatz ([6]) works independent from 
the ratio of healing length and radial size of the BEC cigar. In the limit that the healing length 
is much less than the radial size, TF wave functions are used, in the opposite limit, a Gaussian 
ansatz for the radial part of the wave function 4'n( r ) is appropriate. 

For axial excitations characterized by a wavelength A = 2ir/k exceeding the radial size R± of 
the condensate, we have kR± <C 1, and the dispersion relation reads, in the TF limit for the radial 
wave function [46| 

el k = 2u 2 ± n(n + l) + c 2 k 2 , (7) 

where cq = W /jl/2 and the chemical potential fi = oj\R\/2. Observe that the central speed of 
sound Co is reduced by a factor y/2 from the well-established value ^gp m = ^J~^ for an infinitely 
extended liquid (46| . 

The Eqs. © and ([7J) can be generalized for an expanding condensate. Substituting in Eq. © the 
rescaled radial wavefunction <p n = 4> n (r/b±), and inserting the result into the action dU), integrating 
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over the radial coordinates, we find the following effective action for the axial modes of a given 
radial quantum number n: 



S n = dtdz 



(d tXn ~ v z d z f + c 2 n (z)(d zX n) 2 + M 2 n {z)xl , 



2g 

where the common "conformal" factor C n (z) is given by 



b 2 ± C n (z) = / d'rti. (9) 



r<r„ 



The integration limits are fixed by the z dependent radial size of the cigar = — z 2 /R^b 2 ). 

The averaged speed of sound reads 

= 77T2- / d 2 rp (f) 2 n , (10) 



CrM 



r<r„ 



and the (space and time dependent) effective mass term is, for a given radial mode, obtained to be 

M n(^t) = 7 ^r[ d 2 rp [d r <p n ] 2 . (11) 

^nU± Jr<r m 

The phonon branch of the excitations corresponds to the gapless n = part of the spectrum in 
Eq. ([7]). In this case the radial wavefunction cj)Q does not depend on the radial variable r j4fjj, and 
the mass term vanishes, Mq = 0. We then obtain the following expressions, 

Co(z) = nRl (l - y^j , (12) 

and for the z dependent speed of sound (c = c n= o): 

^ = A( 1_ ^f)- (13) 

We will see below that we need these expressions in the limit z — > only, because only in this limit 
we get the desired exact mapping of the phonon field to a quantum field propagating in a 1+1D 
curved space-time. 



III. THE 1+1D DE SITTER METRIC IN THE CONDENSATE CENTER 

We identify the action ([8]) with the action of a minimally coupled scalar field in 1+1D, according 
to Eq. ([T|). Such an identification is possible only close to the center z = 0, as will now show. The 
contravariant 1+1D metric may generally be written as [l2| 



where A c = A c {x^) is a conformal factor and c = c(z = 0). Inverting this expression, we get a 
covariant metric of the Painleve-Gullstrand form 14711. 



The term y 7 — gg^ contained in the action ([I]) gives the familiar conformal invariance in a 1+lD 
space-time, i.e., the conformal factor A c drops out from the action and thus does not influence the 
(classical) equations of motion. We therefore leave out A c in the formulae to follow, but it needs 
to be borne in mind that the metric elements are defined always up to the factor A c . 

The actions (pQ) and ([8]) can be made consistent if we renormalize the phase field according to 
<3? = Z<3? and require that 

9 c 

holds. The factor Z does not influence the equation of motion, but does crucially influence the 
strength of detector response (see section HVAI below). In other words, it renormalizes the coupling 
of our "relativistic field" $ to the detector. More explicitly, Eq. (fl~6j) leads to 



^=^V^f(^) SB = co, 1S , (17) 

According to the above relation, we have to impose that the expansion of the cigar in the perpen- 
dicular direction proceeds like the square root of the expansion in the axial direction. The constant 
quantity B can be fixed externally (by the experimentalist), choosing the expansion of the cloud 
appropriately by adjusting the time dependence of the trapping frequencies oj\\(t) and uj±(t), as 
prescribed by the scaling equations (i"" 



42] 



Ujf, CO?, w 2 ^2 

gfe- ah- "^^-gd-^n- (18) 



Since both Co and c depend on z, Eqs. (fl2|) and (|T3|) . an effective space-time metric for the 
axial phonons can be obtained only close to the center of the cigar-shaped condensate cloud. This 
is related to our averaging over the physical perpendicular direction, and does not arise if the 
excitations are considered in a higher-dimensional situation, where this identification is possible 
globally. The reason for the 1+lD geometry is of practical origin: It implies that during expansion, 
the coupling g can be chosen independent of time to reproduce the de Sitter metric, which is in 
stark contrast to the 2+1D and 3+1D case, where the interaction needs to increase exponentially 
for a de Sitter effective space-time to be obtained [2(|. Also note that the action (jHJ) does not 
contain a curvature scalar contribution of the form oc XnV~ 6 RlS/wO^)]) i- e -' that it only possesses 
trivial conformal invariance 

We now impose, in addition, the requirement that the metric is identical to that of a 1+lD 
universe, with a metric of the form of the de Sitter metric in 3+ ID (29I . I30I ], We first apply the 
transformation c^dt = c(t)dt to the line element defined by (I15p . connecting the laboratory time 
t to the time variable t. Defining v z /c = vAz = (Bb/co)z (note that the "dot" on b and other 
quantities always refers to laboratory time t), this results, up to the conformal factor A c , in the 
line element 

ds 2 = -Cq(1 - Az 2 )dP - 2c z\/ r Adtdz + dz 2 . (19) 
We then apply a second transformation c^dr = c^dt + z\fh.dzj{\ — Az 2 ), with a constant A. We 



are thus led to the 1+lD de Sitter metric in the form [2£ 

ds 2 = -cl (1 - Az 2 ) dr 2 + (l- Az 2 ) _1 dz 2 . (20) 

The transformation between t and the de Sitter time r (on a constant z detector, such that dt = dt), 
is given by 

f = exp[i?6r], (21) 
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where the unit of lab time to ~ w m 1 is set by the initial conditions for the scaling variables b and b±. 
The temperature associated with the effective metric (|20p is the Gibbons-Hawking temperature 

Hi 

The "surface gravity" on the horizon has the value an = Cq-\/A. = coBb, and the stationary 
horizon(s) are located at the constant values of the z coordinate 




±z H = ±Rn\ (23) 



Combining (f2"2"j) and (f2"3"j) . we see that zu/R\\ is small if uj\\/Tds *C 47r. Therefore, the de Sitter 
temperature needs to be at least of the order of for the horizon location(s) to be well inside 
the cloud. The latter condition then justifies neglecting the z dependence in Co and c in Eq. ([IB]) . 
Though there is no metric "behind" the horizon, i.e. at large z, this should not affect the low-energy 
behavior of the quantum vacuum "outside" the horizon, i.e. in the center of the cloud. 



IV. DETERMINING THE QUASIPARTICLE CONTENT OF THE QUANTUM FIELD 



As stated, the particle content of a quantum field state depends on the observer [1H3J. To 
verify this statement, we need to make the notion of observer operationally precise in that we 
find a quasiparticle detector which couples to the quantum field in question and hence measures 
only the quasiparticles corresponding to excitations above the appropriate vacuum. Specifically, to 
detect the Gibbons-Hawking effect in de Sitter space, one has to set up a detector which measures 
frequencies in units of the inverse de Sitter time r, rather than in units of the inverse laboratory 
time t. The de Sitter time interval dr = dt/bB = dt/y/bb± can be measured by an "Atomic 



Quantum Dot" (AQD) [3ll. |32|. |49||. The measured quanta can then, and only then, be accurately 
interpreted to be quasiparticles coming from a Gibbons-Hawking type process with a constant de 
Sitter temperature (|22p . In addition, we use the tunability to other time intervals feasible with 
our detector scheme, and contrast the thermal de Sitter result with what the detector "sees" if 
tuned to laboratory and "adiabatic" time intervals. We will see that such a tuning is possible in 
the quantum-optical context by temporal variations of the laser amplitudes coupling the hyperfine 
levels of the atom. 

The AQD can be manufactured in an ultracold gas of atoms possessing two hyperfine ground 
states a and (3. The atoms in the state a represent the superfluid cigar, and are used to model 
the expanding universe background on which the phononic quasiparticles propagate. The AQD 
itself is formed by trapping atoms in the state j3 in a tightly confining optical potential V op t. 
The interaction of atoms in the two internal levels is described by a set of coupling parameters 
g c d = 4:ira C( i (c,d = {a,/3}), where a c d are the s-wave scattering lengths characterizing short-range 
intra- and inter-species collisions; g aa = g, a aa = a s . The on-site repulsion between the atoms 
(3 in the dot is U ~ gpp/P, where I is the characteristic size of the ground state wavefunction 
of atoms /3 localized in V op t- In the following, we consider the collisional blockade limit of large 
U > 0, where only one atom of type /3 can be trapped in the dot. This assumes that the double- 
occupancy gap U is larger than all other relevant frequency scales in the dynamics of both the 
AQD and the expanding superfluid. As a result, the collective coordinate of the AQD is modeled 
by a pseudo-spin- 1/2, with spin- up/spin-down state corresponding to occupation by a single/no 
atom in hyperfine state j3. 
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We first describe the AQD response to the condensate fluctuations in the Lagrangian formalism, 
most familiar in a field theoretical context. The detector Lagrangian takes the form 



( -^n* 



nVMMjP ex P 



-i I gp o (0,t')dt' + i5<t> 
'o 



r) + exp 



* / ffPo(O) t )dt — i5(j) 
o 



n). (24) 



Here, A is the detuning of the laser light from resonance, po(z = 0,t) is the central mean-field 
part of the bath density, and I is the size of the AQD ground state wave function. The detector 
variable rj is an anticommuting Grassmann variable representing the pseudo-spin degree of freedom 
of the AQD. The second and third lines represent the coupling of the AQD to the surrounding 
superfluid, where 5(f> and 5p are the fluctuating parts of the condensate phase and density at z = 0, 
respectively. The laser intensity and the effective transition matrix element combine into the Rabi 
frequency below we will make use of the fact that f2 can easily experimentally be changed as a 
function of laboratory time t, by changing the laser intensity with t. 
To simplify (|24j) . we use the canonical transformation 



t] — y Tj exp 



t 

\ I gp (0,t')dt' + it 
'o 



(25) 



The above transformation amounts to absorbing the superfluid's chemical potential and the fluc- 
tuating phase 5(p into the wave function of the AQD, and does not change the occupation numbers 
of the two AQD states. The transformation (|25p gives the detector Lagrangian the form 



d 



L A QD=i( J t V* ) fj-ny/po(0,t)P (fj + fj*) 



-A + (g a p - g)p (0, t) + g af3 5p + —5(f> 



The laser coupling (second term in the first line) scales as 6 _1 / 2 6j 1 , and hence like the de Sitter 
time interval in units of the laboratory time interval, drjdt. We suggest to operate the detector at 
the time dependent detuning A(t) = {g a p — g)po(0,t) = (g a p — g)p rn /(b 2 B 2 t 2 ), which then leads 
to a vanishing of the first two terms in the square brackets of (|26p . 

We now introduce back the wave function of the AQD stemming from a Hamiltonian formula- 
tion, ifi = 4 , p\f3) + Tp a \a). An "effective Rabi frequency" may be defined to be loo = 2£l^J p m l 3 ; at 
the detuning compensated point, we then obtain a simple set of coupled equations for the AQD 
amplitudes 



. dip a _ 
dr ~ 2 



(26) 



where r is the de Sitter time. 

We have thus shown that the detector Eqs. (|26h are natural evolution equations in de Sitter 
time r, if the Rabi frequency is chosen to be a constant, independent of laboratory time t. We 
will see in sections IIVBI and IIV CI that, adjusting in a certain time dependent manner, within 
the same detector scheme, we can reproduce time intervals associated to various other effective 
space-times. 

The coupling of the AQD to fluctuations in the superfluid is described by the perturbation 
potential 



SV(t) = (g a p - g) Bb(r)5p(r). 



(27) 



Neglecting the fluctuations in the superfluid, the level separation implied by (|26p is ujq, and the 
eigenfunctions of the dressed two level system are |±) = (\a) ± \ j3))/^/2. The quantity ujq therefore 
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plays the role of a frequency standard of the detector. By adjusting the value of the laser intensity, 
one can change ujq, and therefore probe the response of the detector for various phonon frequencies. 
Note that if g a p is very close to g, to obtain the correct perturbation potential, higher order terms 
in the density fluctuations have to be taken into account in the Rabi term of (1261) . 

To describe the detector response, we first have to solve the equations of motion (J8|) for the 
phase fluctuations, and then evaluate the conjugate density fluctuations. The equation of motion 
5Sq/5xo = is, for time independent B, given by 

B 2 b 2 d t (b 2 d tX0 ) - 777-^*, {c 2 (zb)Co(z b )d ZbX o) = 0, (28) 

where Z\, = z/b is the scaling coordinate. Apart from the factor Co(z&), stemming from averaging 
over the perpendicular direction, this equation corresponds to the hydrodynamic equation of phase 



fluctuations in inhomogeneous superfluids |50j]. At t — > — oo, the condensate is in equilibrium and 
the quantum vacuum phase fluctuations close to the center of the condensate can be written for a 
given wavevector k in the following form 



Xo = J . n a k exp [-ie 0) fct + ikz] + H.c. (t ->■ -oo), (29) 

y 4Co(UJ/t||eo,fc 

where at are the annihilation and creation operators of a phonon. The intial quantum state of 
phonons is the ground state of the superfluid and is annihilated by the operators fife. With these 
initial conditions, the solution of (1281) is 



l/ 4CoW |eo / fceXP 



■*/ ~B^ +lkZ \ 



+ H.c. (30) 



The solution for the density fluctuations, which are in a superfluid canonically conjugate to the 
phase fluctuations, 5p = —dtxo/g and [Sp(r), Xoir')] = ~~ r ')/{Bb 2 ) [Elj], therefore is 



^ = V^o(om |e o,J (afceXP 



ll Z Bb¥ +lkZb 



+ H.c. (31) 



The Eqs. (|30|) and (|3ip completely characterize the evolution of the quantized hydrodynamic con- 
densate fluctuations when the transverse modes are not excited. Observe that the evolution pro- 
ceeds without frequency mixing in the adiabatic time interval defined by dr a = dt/Bb 2 (the "scaling 
time" interval dt/B 2 b 2 defined in [2(| is proportional to this adiabatic time interval). Therefore, 
in the "adiabatic basis," no frequency mixing occurs and thus no quasiparticle excitations are cre- 
ated, see also section HV CI below. This hints at a hidden (low energy) symmetry, in analogy to the 



(exact) 2+1D Lorentz group SO(2,l) for an isotropically expanding BEC disk, discussed in |44|]. 



A. Quasiparticle detection in de Sitter time 

The coupling operator 5V causes transitions between the dressed detector states |+) and |— ) 
and thus can be used to effectively measure the quantum state of the phonons. We consider the 
detector response to fluctuations of ^, by going beyond mean-field and using a perturbation theory 
in 5V . There are two physically different situations. The detector is either at t = in its ground 
state, (|a) + \P))/y/2, or in its excited state, (|a) — \f3))/^/2. We define P+ and P_ to be the 
probabilities that at late times t the detector is excited respectively de-excited. Using second order 
perturbation theory in SV, we find that the transition probabilities for the detector may be written 
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where the absolute square of the transition matrix element is given by 



\ T ±\ = 

Calculating the integrals, we obtain 

P± = J 



dr 

W) 



exp 



±ieo,k 



dr' 



77 + iuqt 
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2BbR\\Co{0) I 1 + n B 
where a (formally divergent) sum is contained in the factor 

k eo > k 



(33) 



(34) 



(35) 



and we obtain Bose-Einstein distribution functions at boson energy ujq, containing the de Sitter 
temperature (f22|) : 



riB 



exp[w /r dS ] 



(36) 



We conclude that an expansion of the condensate in z direction, with a constant rate faster than 
the harmonic trap oscillation frequency in that direction, gives an effective de Sitter space-time, 
characterized by a thermal distribution at the de Sitter temperature T^s, as measured by the AQD 
in its natural time interval dr. 

We now show that J in f)35[) is proportional to the total de Sitter time of observation, so that 
the probability per unit time is a finite quantity [52j]. At late times, the detector measures phonon 
quanta coming, relative to its space-time perspective, from close to the horizon, at a distance 
Sz = zh — z <C zh = A" 1 / 2 . The trajectory of such a phonon in the coordinates of the de Sitter 
metric ([20]) . at late times r, is given by 



In 



£H 

5z 



2\/Acot. 



(37) 



This implies that the central AQD detector measures quanta which originated at the horizon with 
the large blueshifted frequency 



e 0,k 



—^L exp[coT\/A]. 



(38) 



Making use of the above equation, we rewrite the summation over k in (|35p as an integral over 
detector time: 



eo,fc 



Ruljq f de 



7TC 

vrc 



-O.Jfc 



RnLJo 

7TC 



VXco J dr 



Bbr . 



(39) 



Therefore, the probabilities per unit detector time (de Sitter time) read, where upper/lower entries 
refer to P + /P_, respectively: 



dP± 



9a/3 



1 B 



2 9^o 



2C (0)c 1 l + n B 



(40) 
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They are finite quantities in the limit that r — > oo. In laboratory time, the transition probabilities 
evolve according to 



P ± (t) = P ^ln 



, (41) 
1 + n B ' V ' 



where, from relation (flfi]) . Pn = Z 2 [(g af} /g - l)B} 2 /2. We see that the detector response is, as it 
should be, proportional to Z 2 , the square of the renormalization factor of the phase fluctuation 
field. 

The absorption and emission coefficients dP±/dr satisfy Einsteinian relations. Therefore, the 
detector approaches thermal equilibrium at a temperature Tas on a time scale proportional to 
Z^Uq 1 . Remarkably, our de Sitter AQD detector thus measures a stationary thermal spectrum, 
even though its condensed matter background, with laboratory time t, is in a highly nonstationary 
motional state. Since Z 2 oc \J p m o? s (uj±/p) 2 , not-too-dilute condensates with uj± ~ /i (i.e., close to 
the quasi-lD regime [H^]) are most suitable to observe the Gibbons-Hawking effect. 

The verification of the fact that a thermal detector state has been established proceeds by 
the fact that the two hyperfme states a and (3 are spectroscopically different states of the same 
atom, easily detectable by modern quantum optical technology. When the optical potential is 
switched on, the atoms are in the empty a state originally, which is an equal-weight superposition 
of |+) = (|a) + |/3))/V2 and |-) = (\a)-\/3))/y/2. The thermalization due to the Gibbons-Hawking 
effect takes place in the dressed state basis consisting of the two detector states, i.e. of the states 
|+) and |— ), on a time scale given by the quantities P± in Eq. (|4ip . For the laboratory observer, 
the Gibbons-Hawking thermal state will thus appear to cause damping of the Rabi oscillations 
on the thermalization timescale, i.e., friction on the coherent oscillating motion between the two 
detector states occurs, due to the thermal phonon bath perceived by the detector. The occupation 
of the detector states can be measured directly using atomic interferometry: A 7r/2-pulse brings 
one of them into the filled {(3) and the other into the empty (a) state. To increase the signal 
to noise ratio, one could conceive of manufacturing a small array of AQDs in a sufficiently large 
cigar-shaped host superfluid, and monitor the total population of (3 atoms in this array. 



B. Detection of quasiparticles in laboratory time 

We contrast the above calculation with the response the AQD detector would see if tuned 
to laboratory time. This can be realized if we let f! oc i, such that the quantity fiy po{0, t) = 
Qy/Pm/ (Bbt) = const., in the Rabi term on the right-hand side of (|26p . is time independent. 
The detector has, therefore, the laboratory time interval dt in the equations for the occupation 
amplitude with this choice for Q = Q(t). The Painleve-Gullstrand metric (|15|) in pure laboratory 
frame variables, assuming B 2 b 3> u 2 /b 2 like in the derivation of the de Sitter metric (|20p . reads 

ds 2 = %— (l - Az 2 ) dt 2 - —dzdt + dz 2 . (42) 

B 2 b 2 t 2 y J t v ' 

The metric (|42p is asymptotically, for large t, becoming that of Galilei- invariant ordinary Newtonian 
laboratory space, i.e. is just measuring length along the z direction, because the speed of sound in 
the ever more dilute gas decreases like 1/t and the "phonon ether" becomes increasingly less stiff. 
The transition probabilities for absorption respectively emission are now given by 
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where the matrix elements are, cf. Eq. ([3 

dt 



Bb 2 



cxp 



dt' 



+ iuj^t 



Substituting the adiabatic time interval dr a 



dt/{Bb 2 t 2 ) leads for large t to 

■ - (BbHy 1 , 



where tq & = dt(Bb 2 t 2 ) 1 . The transformation to adiabatic time maps t G 
r a € [— oo,Toa] and, by further substituting y = eo ) fc(r a — ro a ), we have 

2 



\T±\ 



1 



' OA' 



dyexp 



o 



y=F 



^QfOfc 1 



(44) 
(45) 

-oo, +oo] onto 
(46) 



The integral is a linear combination of Bessel functions. To test its convergence properties, we 
are specifically interested in the large eofc limit. Performing a stationary phase approximation 
for large A = iwoeofc/ Bb 2 , we have for positive A (absorption) that the integral above becomes 
J (A) = (7r\/^4) 1//2 exp[-2\/]4] and for negative A (emission) J (A) = (7r^/\A\) 1/2 . The final result 
for the laboratory time transition probabilities then is 



p± 




Bb 2 



1 



9a/S 

(47) 

where Ep\ ~ /i is the ultraviolet cutoff in the integral for the emission probability P- , the "Planck" 
scale of the superfluid. Because of the convergence of the absorption integral for P + , the total 
number of quasiparticles detected remains finite, and there are thus no quasiparticles detected by 
the laboratory-frame-tuned detector at late times. This is in contrast to the de Sitter detector, 
which according to (jlU]) detects quasiparticles in a stationary thermal state. 



C. Adiabatic "Minkowski" basis with no quasiparticles detected 

There is a detector setting which corresponds to a detector at rest in the Minkowski vacuum. 
This setting is represented by the adiabatic basis, with time interval defined by dr a = dt/Bb 2 , 
realizable with the AQD by decreasing the Rabi frequency linearly in time, oc 1/t. Then, no 
quasiparticles whatsoever are detected at any instant, i.e., no frequency mixing of the positive and 
negative frequency parts of (f3T|) does take place. The associated space-time interval 

ds 2 = b 2 {-cldrl + dzl\ (48) 

is simply that of (conformally) flat Minkowski space in the spatial scaling coordinate and adi- 
abatic time coordinate r a . Specifically note here that the adiabatic time interval is not equal to 
laboratory time, (fr a oc dt/t 2 . The detector in its effective Minkowski space hence measures, us- 
ing this proper time interval, that dissipation-free expansion of the BEC cigar takes place. The 
"absolute" time of the experimentalist in the lab frame can thus be different from the proper time 
interval of the detector in the sonic Minkowski frame. The "internal" observer therefore can as- 
certain independently from the lab experimentalist that indeed dissipation-free expansion of the 
background takes place and no real quasiparticles are created. 
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V. CONCLUDING CONSIDERATIONS 



We have demonstrated that the dependence of quasiparticle detection on the motional state of 
the detector and on the way the detector couples to the quantum field can be made experimentally 
manifest in a tabletop setup. An atomic quantum dot placed at the center of a linearly expand- 
ing cigar-shaped Bose-Einstein condensate has a proper time interval which can be tuned using 
quantum optical means like the intensity of a laser. The proposed quantum optics experiment can 
thus confirm in the effective curved space-time setting of analogue gravity that indeed "A particle 
detector will react to states which have positive frequency with respect to the detector's proper 
time, not with respect to any universal time []J|." 

The equilibration time scale of the detector, and thus the time scale on which the Rabi oscilla- 
tions between the detector states are damped out, is set by the detector frequency standard (the 
level spacing) cjq, and by the renormalization factor Z: 



?~equil 



Z- 2 ^ 1 oc (p m al)~ 1/2 (/Vw±) 2 ^o" 1 - (49) 



The renormalization factor Z contained in (I17D determines the equilibration rapidity because it 
physically expresses the strength of detector-field coupling, resulting in the transition probabilities 
(|41j) . It is related to the initial diluteness parameter D p (0) = (^dj) 1 / 2 of the Bose-Einstein 
condensate and to the ratio fJ,/u)j_, which determines inasmuch the system is effectively one- 
dimensional, Z 2 oc D p (0)(uj±/ fi) 2 . To obtain sufficiently fast equilibration, the condensate thus 
has to be initially not too dilute as well as close to the quasi-lD regime, for which the trans- 
verse harmonic oscillator energy scale is of order the energy per particle, [x ~ lo±. These two 
conditions have another important implication. The ratio of the instantaneous coherence length 
£ c (t) = (871730(0, t)a s )~ l l 2 oc t and the location of the horizons z = zu = ±A -1 / 2 , which are sta- 
tionary in the present setup, has to remain less than unity within the equilibration time scale. If 
this is not the case, the coherence length, which plays the role of the analogue Planck scale (which 
is laboratory time dependent here), exceeds the length scale of the horizon at equilibration, and 
the concept of "relativistic" phonons propagating on a fixed curved space-time background with 
local Lorentz symmetry becomes invalid. The ratio £ c (t)/zH at the lab equilibration time scale 
t = tequil = *o exp[27r(Tds/^o) (/V^J-) 2 Dp 1 (0)] following from the de Sitter equilibration time in 
Eq. (I49p . expressed in experimentally relevant scaled parameters, reads 



£c(*cquil) ^O^dS 

— = — exp 

m p m a s 



2vr^ (-!±X 1 
wo \u±J D p (0) 



(50) 



We see that this ratio changes exponentially with both the initial diluteness parameter D p (0) 
and the quasi-lD parameter fi/u;±. In most currently realized Bose-Einstein condensates, with 
comparatively long lifetime, the diluteness parameter D p ~ 10~ 2 . Here, we initially need D p (0) < 
0(1) to have the condition £ c (i C quii) / zh < 1 fulfilled, assuming a reasonably large value of the 
de Sitter temperature T^s- Though the condensate has to be initially quite dense, it is to be 
stressed that the central density decays like t~ 2 during expansion. Therefore, the rate of three- 
body recombination losses quickly decreases during the expansion of the gas, and the initially 
relatively dense Bose-Einstein condensate, which would rapidly decay if left with a D p close to 
unity, can live sufficiently long, the total rate of three-body losses decreasing like po(0,i) oc t~ 4 . 
The gas parameter D p plays the roleof a loop expansion parameter, which measures the extent 



to which a mean-field theory applies [4l|]. The Gibbons-Hawking effect in the BEC is therefore 
intrinsically quantum because the signal measured by the detector, Eq. (|4T|) . as already stated 
above, contains this "dimensionless Planck constant" D p (0) at the initial laboratory time. When 
D p becomes large, of order unity, the mean-field theory employed in the above considerations 
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strictly speaking breaks down. Because the gas parameter enters the magnitude of the signal, a 
sufficiently large signal-to-noise ratio, together with rapid equilibration of the detector, can however 
be achieved only by using initially dense clouds with strong interparticle interactions and relatively 
large gaseous parameter. The required initial state can be prepared by using a sufficiently rapid 
external magnetic field sweep to a Feshbach resonance to increase a s for an already existing BEC, 
and shortly thereafter (to avoid large three-body losses) expand, according to the required temporal 
behavior of scaling parameters in Eq. (I18p . 

Considering the nonclassicality aspect of the Gibbons-Hawking effect from a different angle, 
the phononic quasiparticles of the superfluid can be regarded as non-interacting only in a first 
approximation in D p <C 1. Self-interaction between the phonons, induced by larger values of the 
gas parameter, can lead to decoherence and the relaxation of the phonon subsystem. The same 
line of reasoning applies to the evolution of interacting quantum fields in the expanding universe. 
The interactions between quasiparticle excitations in the BEC induced by them scattering off each 
other can therefore be related to decoherence processes in cosmological models of quantum field 
propagation and particle production. Finally, as already mentioned in section m the Einsteinian 
part of the effective action becomes increasingly important for larger D p , so that the combined 
effects of backaction of quantum fluctuations on the sonic metric and of decoherence due to phonon- 
phonon interactions may lead to increasingly rich behavior of the quantum field evolution. 
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